Google 



This is a digital copy of a book that was preserved for generations on library shelves before it was carefully scanned by Google as part of a project 

to make the world's books discoverable online. 

It has survived long enough for the copyright to expire and the book to enter the public domain. A public domain book is one that was never subject 

to copyright or whose legal copyright term has expired. Whether a book is in the public domain may vary country to country. Public domain books 

are our gateways to the past, representing a wealth of history, culture and knowledge that's often difficult to discover. 

Marks, notations and other maiginalia present in the original volume will appear in this file - a reminder of this book's long journey from the 

publisher to a library and finally to you. 

Usage guidelines 

Google is proud to partner with libraries to digitize public domain materials and make them widely accessible. Public domain books belong to the 
public and we are merely their custodians. Nevertheless, this work is expensive, so in order to keep providing tliis resource, we liave taken steps to 
prevent abuse by commercial parties, including placing technical restrictions on automated querying. 
We also ask that you: 

+ Make non-commercial use of the files We designed Google Book Search for use by individuals, and we request that you use these files for 
personal, non-commercial purposes. 

+ Refrain fivm automated querying Do not send automated queries of any sort to Google's system: If you are conducting research on machine 
translation, optical character recognition or other areas where access to a large amount of text is helpful, please contact us. We encourage the 
use of public domain materials for these purposes and may be able to help. 

+ Maintain attributionTht GoogXt "watermark" you see on each file is essential for in forming people about this project and helping them find 
additional materials through Google Book Search. Please do not remove it. 

+ Keep it legal Whatever your use, remember that you are responsible for ensuring that what you are doing is legal. Do not assume that just 
because we believe a book is in the public domain for users in the United States, that the work is also in the public domain for users in other 
countries. Whether a book is still in copyright varies from country to country, and we can't offer guidance on whether any specific use of 
any specific book is allowed. Please do not assume that a book's appearance in Google Book Search means it can be used in any manner 
anywhere in the world. Copyright infringement liabili^ can be quite severe. 

About Google Book Search 

Google's mission is to organize the world's information and to make it universally accessible and useful. Google Book Search helps readers 
discover the world's books while helping authors and publishers reach new audiences. You can search through the full text of this book on the web 

at |http: //books .google .com/I 



ON PARTIAL DIFFERENTIAL EQUATIONS. 



PRINTKD BY a J. OLAT. MJL 
AT THB UKIVBBSITT FKEUS. 




ELEMENTAEY TREATISE 

ON 

PARTIAL DIFFERENTIAL EQUATIONS 

Ststtjpirt fat ttit 8X)U of jKtulinrttf tn tjtit SSnibmfts. 
GEORaE BIDDELL AIEY, M.A. LL.D., D.C.L. 




lonbon antr tJCamtidtigt : 

MACMILLAN AND CO. 

1866. 

/6JL. J^. ^ 



PREFACE. 



The work now offered to the University is strictly an 
Elementary Treatise. No attempt has been made to go 
into all the varied details, of methods and examples, 
which present themselves in the wide field of Partial 
Differential Equations, considered purely as an Algebraical 
subject 

I have endeavoured, however, to omit no important 
consideration affecting the Principles of those Equations. 
And I trust that the methods of solution here explained, 
and the instances exhibited, will be found sufficient for 
application to nearly all the important problems of Physical 
Science, which require for their complete investigation the 
aid of Partial Differential Equations. 

G. B. AIRY. 



BoTAL Obsebvatobt, Gbsskwioh, 
J 866, Augtut 15. 
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ON PARTIAL DIFFERENTIAL 

EQUATIONS. 



PRELIMINARY NOTICE ON INTEGRATION. 

1. In all that follows, we shall suppose that it is 
always possible to effect simple integration; inasmuch as 
any difficulties of integration, connected with the solutions 
of partial differential equations, do not affect the principle 
of those solutions. Thus, for instance, we shall not 
hesitate to represent an unknown function of x by x' (^)j 
(the second differential-coefficient of %(ic)), on the assump- 
tion that, whatever be the form of x' {x), we can in some 
way find the function 'x^{x) of which it is the second differ- 
ential-coefficient. 

CHARACTERISTICS OF THE SOLUTIONS OF SIMPLE 
DIFFERENTIAL EQUATIONS. 

2. Before entering on the subject of Partial Differ- 
ential Equations, it may be convenient to consider some of 
the characteristics of the solutions of Simple Differential 
Equations. 

3. To begin with Simple Differential Equations of 
the first order. Suppose, for facility of geometrical illus- 

^ ^. B 
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tration, we consider the equation y-^^o,^ or the equation 

y "j- ^^y <>f which the former, translated geometrically, 

indicates that the subnormal of a plane-curve (to be found) 
is constant, and the latter indicates that the subtangent is 
constant. The algebraical solutions are easily found: in 
each, there is a constant, which does not occur in the 
original differential equation, and is not defined by it; a 
constant of that class described (perhaps improperly) by 
the term " arbitrary," but which really means " not yet 
determined, but enabling us by proper determination of its 
value so to fix the value of a: corresponding to a given value of 
y that we can adjust the solution to some specific condition." 
The reader is requested to observe that instead of the 
term "arbitrary constant, "we shall always use the term 
"undetermined constant." 

4. If we treat the geometrical translation of the 
differential equation by a geometrical process, always 
drawing a normal or tangent so as to make the subnormal 
or subtanpnt constant, then drawing by means of it a 
small portion of the curve, then repeating the process, &c., 
we may produce a polygon which will approach to the 
strict solution, with smaller errors (by taking the sides 
small enough) than any small quantity that can be assigned. 
In each case, however, a starting point is necessary. 

5. In both ways of treating the problems these con- 
ditions manifestly hold : — 
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The curve, in each case, is one definite curve. 

The curve, in each case, is a continuous curve, ex- 
pressed by the same equation through its whole 
extent. 

Even if there be isolated points or curves, still the 
same one equation defines the whole. 

It is necessary to introduce one undetermined quantity, 
enabling us to adjust the curve to a specific con- 
dition given by special considerations: that unde- 
termined quantity is however a simple constant. 

6. Let us now consider a Simple Difierential Equation 
of the second order: such for instance as is given by this 
problem, "To find the curve in which the radius of curva- 
ture is a function, given in form, of the ordinate." Here 
the algebraic solution gives a formula requiring two unde- 
termined quantities, still simple constants: the equivalent 
geometrical treatment shews that we require two elements, 
(as for instance, the value of x and the inclination of the 
tangent, for some one value of y). But all the conditions 
hold which are mentioned in Article 5 : the only difierence 
being that, instead of one undetermined quantity, a 
simple constant, there must now be two undetermined 
quantities, simple constants. 

7. The conditions which will be found to hold in the 
solutions of Partial Differential Equations differ very 
remarkably from these* 



B2 



4 OH PABTIAL DIFFERENTIAL EQUATIONS. 

IKTBODUCnON OP TWO OB MORE IXDEPEXDENT 

VARIABLES. 

8. It in convenient, in Simple DiflFerential Equations, 
to contrjmplate one quantity z as expressed algebraicallj 
in terms of another quantity x (whether it be actually so 
expressed, or not): and to consider its differential-coefficients 

dz d fdz\ d^z p 
dx ' dx\dx) da? ' ' 

as b(5ing formed by taking the limiting values of fractions, 
in each of which the denominator is an increment in the 
value of a?, and the numerator is the corresponding incre- 
ment in tlie last preceding term of the series. And this is 
cxprosflod by saying that x is the "independent variable." 
(ieomotrically 'it is illustrated by supposing that we consider 
on ordinate « of a plane curve, as also the inclination-tan- 
gent of the curve's-tangent, the change of that inclination- 
tangent for a small change in the value of a?, &c., to be 
expressed in terms of the abscissa x. 

0. There is no difficulty in conceiving z to depend on 

two quantities x and y, combined in any way and with any 

constants under any functional formula; and in considering 

that wo may at pleasure change the value of x without 

changing that of y, or may change the value of y without 

changing that of a?, or may change both simultaneotv^ly. 

If wo do not vary y, y is pro tempore a constant: and by 

dz d^z 
varying x alone, we may form -r- > -i-t > &c., just as in the 



ILLUSTRATION FEOM SOLID GEOMETRT. 5 

functions which we have illustrated by reference to a plane 
curve. If we do not vary a?, then x is pro tempore a con- 

(iz d z 

stant, and by varying y alone, we may form -^ , -y-g , &c. 

But if we vary both x and y, then we have both series of 

cPz 

differential-coefficients lust set down, and also ^ 7- or 

•^ ax . ay 

-J — -J- (which, as is known in the expansions obtained by 

the Differential Calculus, have the same value), with other 
coefficients of analagous form, which do not appear in such 
series as those of Article 8. Here we consider x and y as 
" two independent variables," which are strictly indepen- 
dent of each other : and 2? as a function of both. 



10. Analytical Geometry of three dimensions assists 
greatly in illustrating these algebraical conceptions. If cc, 
y, «, be three rectangular co-ordinates, then the expression 
of z by means of x and y determines the numerical value 
of the height of an ordinate z which must be erected over 
the point on the plane x, y, defined by the numerical value 
of X and y, in order that the elevation of the summit of that 
ordinate z may represent the value of our algebraical func- 
tion. If we do this for an infinite number of values of x 
and y, we determine an infinite number of ordinate-summits, 
all which lie in one curved-surface. If, supposing this 
curved-surface formed, we then contemplate the values of z 
with X invariable and y variable, we shall include all 
these values which are in a plane parallel to the plane yzj 
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and at the given distance x from the origin of co-ordinates. 

This is the same as forming a section of the curved-surface 

by a plane parallel to the plane yz, and considering the 

curved-intersection as an ordinary plane-curve. Thus we 

dz cPz 
can obtain -^ , ^-s , &c. as in Article 8 : but it must be 
dy' df' 

remembered that, in general, x is lurking in their expres- 
sions as a constant; and therefore, in the complete differenti- 
ation (with respect to x) of any formula in which they may 
occur, each must be considered as a function both of x and 
of y. Similarly, the values of z with y invariable and x 
variable will be formed by a section of the curved-surface, 
made by a plane parallel to ocz : the differential-coefficients 

dz d^z 

-T- , ^-5j, &c. can be formed; and a similar caution with 

regard to their dependence both on x and on y must be 
borne in mind. 

11. The coefficients -r- , -r-^ , &c. are, as in Article 

8, to be considered as defining the inclination-tangent 
of the curve's-tangent, the change of that inclination- 
tangent for a small change in the value of x, &c., &c., in 
the curve formed by the intersection of the plane parallel to 
xz with the curved-surface, lu like manner, the coefficients 

dv ' d7 ' ^^'' define similar elements in the curve formed 

by the intersection of the plane parallel to yz with the 

d^z 

curved-surface. But the coefficient -j j- requires further 

ax . ay 
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explanation. Its real expression is j- [--f] • Now -=- de- 
fines the inclination-'tangent of the curve' s-tangent in the 
sectional plane parallel to yz : and therefore its differential- 
coefficient with respect to x defines the rate at which that 
inclination-tangent, in the plane parallel to yz, is changed 
Ly shifting the sectional-plane in the direction of x. 
Analogous explanations apply to analogous succeeding 
terms. 

12. Algebraically, there is no difficulty in conceiving 
£ as a function of any number of independent variables, 
as Uy t7, Wf X, y, and in treating the various differential co- 
efficients to any degree of complexity. But we cannot 
extend any further the elucidation derived from Solid 
Geometry. 

13. Although, in strictness, no difference is really 
made by any physical peculiarity in the nature of the 
quantities represented by the independent variables, pro- 
vided that their magnitudes can be defined by numbers and 
can therefore be treated by algebraical process ; yet it may 
be well to mention that, in some of the most important 
applications of Partial Differentials, one of the independent 
variables is the expression for time. Thus, suppose that 
we are considering the nature of the disturbance which 
each particle of air is undergoing in a musical pipe: we 
want to know the entire movements of the particle whose 
original ordinate was X; and this implies that, in the 
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general formula for the disturbance, x (the general value of 
X) is to be made constant for that particle, but that no 
limit is to be put on the variations of t We also want to 
know the state of disturbance of all the particles at a 
certain time T, and this implies that t is then to be made 
constant, but that x is to admit of all possible variation. 
And these can only be comprehended in a general formula 
which admits of variations both of x and of t Similarly 
for the motion of the particles of waves of water, &c. &c. 

14. As, in Simple Differential Equations, the data of 
the problem, whether mechanical or geometrical, usually 

lead to an equation between a?, ^? ;j~ > 'j~2 > &c., from 

which we desire to obtain a general expression for z in 
terms of a?; so in Partial Differential Equations the data of 

dz 

the problem usually lead to an equation between x,y, z, -j-^ 

dz d^z d\ d^z g ^ , . , J • X 1.x • 

jT-j 3-2> ~7 r > ~7-^> <J^c, from which we desire to obtam 

ay do^' dx . dy djf 

Si general expression for z in terms of x and y. It is the 
object of the present treatise to shew how this may be 
done, and how the meaning of the solutions may be ex- 
plained, in some of the simpler cases : and we now proceed 
with the Solution of Partial Differential Equations. 
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TREATMENT OP THE SIMPLEST PARTIAL DIFFERENTIAL 
EQUATION OP THE FIRST ORDER. 

15. The simplest Partial Differential Equation of the 

first order is, 

dz dz dz dz . 

-^- =a -7-, or -:= — a-7- = 0, 

ax ay ax ay 

[Here and in all subsequent operations we shall put 
uioT ax-\- y.] 

Evidently this equation is satisfied by making 2? = w. 

For, -7- = a, ^- = 1, and therefore -7 a^- = 0\ But 

dx ay ax ay 

it is also satisfied by making z = ^(w), where <f> expresses 
a function whose form is totally unlimited. For then 

^ =*(")-^= ^(«)xa; 

dz ,. . du ,,. V ^ 

5y ='^(")-^ =.^(«)xl; 

and therefore 

^-"^ = f («).{«-«} = 0. 

16. We have thus found a solution of considerable 
generality in form; but there is no evidence that the 
generality is complete. In order to obtain complete and 
certain generality, we shall use a process of Change of 
Independent Variable, guiding ourselves in the selection of 
a new Independent Variable by the indications derived 
from the last imperfect process. 
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17. Instead of considering 2; as a fanction of x and y, 
let us consider 2; as a function of x and u. [Sach a sup- 
position is certainly competent to represent z^ because, if 
we had z expressed by x and y, we have merely to sub- 
stitute for J/ its value u — ax, and z then becomes expressed 
by X and u.] Biaving z then a function of x and u, where 
u is itself a function of x and y, we proceed to express the 
original equation in terms of the new differential-coefficients. 

We shall put -^ for the entire differential-coefficient in 

regard to every, way in which x could appear when x 

dz 
and y were used, and -i- for the differential-coefficient 

when x and u are used: the original equation is therefore 
to be written 

diz^ dz 
ax ay 

18. « is a function of x and m, where u dejpends on 
X and y. 

Therefore 



d{z) ^ dz dz du _ dz dz 
dx dx du' dx dx du ' 



and 



dz ^ dz du ^ dz 
dy du ' dy du 

Substituting these in the equation which is to be solved, 

d(z) dz dz 
dx ay ax 
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* 

That is to say, when z is expressed in terms of x and w, 

dz 

—- = 0. 
dx 

The integration of this simple quantity involves the 
most important considerations of the whole theory. 

19. By the ordinary rules of integration, 2f = (7, where 
G is something which does not contain a?. But what may 
it contain? It may contain every thing whatever except x. 
It may contain any constants whatever. It may contain u\ 
for the differentiation with respect to x will not touch m, 
which is another Independent Variable and is a constant 
with respect to x. It may contain any function of w. It 
may contain these included in the form 

"any function of constants and w." 

But, in ordinary algebraical language, this would be 
sufficiently described as "any function of w," provided that 
we bear in mind that u may even actually disappear, so as 
to leave constants only. This being understood, the solu- 
tion of the equation is, 

z = <^(m.), 

OT z = (j) {ax + y)j 

where (f> expresses a function of any form whatever which 
considerations not yet presented to us may induce us to 
adopt. 

The steps of this process leave no opening for further 
generality in the solution: which is therefore complete. 
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20. Instead of assnming 2; to be expressed hj x and u , 
with elimination of y, we might have assumed 2; to be 
expressed bj y and w, with elimination of x. Then the 
process is, 

dlz) _ dz dz du _ dz dz 
dy dy du ' dy dy du 

dz dz du dz 



dx du' dx du ^ 



and 



r. _ dz d(z) _ dz ^ 
dx dy dy' 



from which the same conclusion follows, namely 

z = (l>(u) = <f>{ax + y). 

21. We may even change the Independent Variables 
so as to eliminate both x and y; and, as this process will be 
important hereafter, we shall shew its application here. 
Let u=^(xx + y, v = ex+fy, where e is not equal to af; 
and let z be supposed to be expressed in terms of u and r. 



Then 



dz _ dz du dz dv _ dz dz 

dx du ' dx dv' dx du dv * 



dz ^ dz du dz dv ^ dz < , ^ /. 
dy " du' dy dv' dy " du dv '' ' 
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and the equation 

dz ^^ _ A 
dx dy ^ 

becomes 

whence, in the same manner as before, 

z = ^(w) = ^ {ax +y). 

GEOMETRICAL INTEBPRETATION OP THE SOLUTION. 

22. In order to shew the geometrical meaning of this, by- 
reference to such a curved-surface as has been considered in 
Articles 10 and 11, we may observe that if Bx and 8y are 
exceedingly small increments of x andy, then the increment 

- dz ^ dz ^ , . , . dz dz , 

01 z = -^ox + -r~oy; which, since ;T- = <^y- » becomes 

dz 

•J- {aSx + Sy) ; and, if we always make Sy = — aSx, the 

tL^jnt of . will U.0: .ha. U, if we alway. ^ fto,„ 
one point on the plane xy to another point on that plane by 
a short line making with the axis of x the angle whose 
trigonometrical-tangent is a (which process continually 
repeated wiU produce a finite straight line in the same di- 
rection), the value of z will be invariable. Therefore, the 
curve-surface will consist of a series of parallel lines, in each 
of which the values of z are the same so long eis ax+y is the 
same, but which is subject to no other condition. It will 
therefore be such as is represented in figure 1 ; a cylindrical 
surface (in the widest sense of the word), whose axis is 
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parallel to xy and makes with x the angle whose trigono- 
metrical tangent is a, but whose transversal section, or 
section by either of the planes xz, yz, is absolutely undefined 
by the equation. Those sections may be adapted to specific 
physical data, and are therefore of the nature of the 
quantities which we call undetermined. But, as the only 
condition established is that the lines (of which the surface 
consists) shall be parallel, the sections need not to consist 
of curves defined each by a single equation in its whole 
extent. A section may consist of a bit of a circle with a 
bit of a cissoid and then a bit of a logarithmic spiral; and 
these may be joined at any angles; or the section may be 
utterly undefinable by algebra: and after it has gone on to 
some extent, there may be an absolute interruption and 
then it may begin again, &c. 

CHARACTEBISTICS OF THE SOLUTIONS OF PARTIAL 
DIFFERENTIAL EQUATIONS. 

23. Thus it appears that these conditions hold in the 
solution ; 

It is not certain that one symbolic expression will 
represent the solution, or that one curve defined by 
one equation will represent the section of the curve- 
surface. 

It is not certain that the section is continuous, either as 
to actual connexion of parts, or as to gradual change 
of direction, or in any other circumstance. 
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There may be isolated parts of the curve, of totally 
distinct character. 

An undetermined quantity must be introduced : but this 
undetermined quantity will not usually be a simple 
constant, but will be (in the most general case) a 
discontinuous function, such as corresponds to the 
discontinuity of curve mentioned in Article 22. 
On comparing these conditions with those in Article 5, 
the reader will see the characteristic diflferences between 
the solutions of Simple and of Partial Diflferential Equa- 
tions. 

INTERPRETATION OP THE SOLUTION WHEN ONE 
INDEPENDENT VARIABLE IS TIME. 

24. There is another instance of the solution in 
Article 19, &c. or (more frequently) of a solution of the 
same form in the treatment of Partial Differential Equa- 
tions of the Second Order, which merits explanation. Let 
one of the Independent Variables be time: put t for ^, and 
suppose the function of t + ax to be '>^ {h {t •\' ax)] or 
^ (J^ + caj), c being = la. Then z^^^^bt+cx). Let any 
one value of ht-\-cx\>Q Bi then if t be increased by c(7 
(where G is really arbitrary), and if a? be increased by 
— J (7, the value of lt-\-cx is still =5, and therefore z 
remains the same. That is, at a time later by c(7, we shall 
have the same value of z provided that we take an abscissa 
smaller by hC This is just the same as if the ordinate z 
was slid backwards upon the abscissa x with a velocity 
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represented hj -. The same applies to every other 

ordinate z. And thus it appears that the whole curve 
defined by the summits of z may be conceived as changing 

its shape by sliding backwards with a definite velocity - . 

c 

This is the characteristic of a wave, in which (as in waves 
of water, to take the most familiar instance) the form 
travels while the particles do not travel with the form. 
Whenever therefore we meet with such a term as (f> {ax+y) 
we may at once take for granted that one of its applica- 
tions may be to the motion of a wave. 

If the value of z had been yjr (bt — cx), we should have 
found in like manner that the wave is travelling forwards. 

These remarks do not in any degree interfere with 
those of Article 23. 



OTHER PARTIAL DIFFERENTIAL EQUATIONS OF THE 

FIRST ORDER. 

dz dz 

26. We proceed with the equation ^- "" ^ t~= ^ (^>y)> 

where the form of the function a is given. We shall solve 
it by the method of Change of Independent Variables. 
Let w = aa; + y, t? = eaj+j^, where e is not equal to af, 
(which, by giving different values to e and /, includes the 
three changes in Articles 19, 20, 21). The value of x in 

terms of u and v is '^ — , that of y is — -^ . Treat- 
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ing the equation as in Article 21, we find 

(e-af)— = a(x v) = a f -^" '"^ Z£?i±^V 
^« «/;^^ - a^.y; - '^\af-e' af-e )' 

whence, integrating on the same principles as before, 

, . 1 [ fju-v -eu+av\ 

where, in the integration with respect to v, u is to be con- 
sidered constant; and, after the integration, ax +y and 
^^ +fy are to be substituted for u and v. 

26. On trying this process upon any ordinary function 
(as a function in integral powers of x and y, or a circular 
or exponential function) of which tlie integrations are easy, 
it will be found that e and/ disappear from the result. It 
may be well to shew this in a simple instance. 

Let one term of a (a?, y) be a;^ . y*. 

• 

The corresponding term of the integral is 

Integrate by parts, beginning with the first factor; it 
gives 

+ 1 



(JP +1) . («/-^) 



_^\p+«+i 



{fu-vy\{-eu+avy 



^3^^/(>-r.(— +-r 



(p+l).(a/-e) 
D. E. 
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and, by repeating the process, we shall obtain the result in 
y +1 terms, in powers of x and y, and without e and/. 

27. If we had begun the integration-bj-parts with 
integration of the second factor, we should have obtained a 
different series of terms of powers of x and y\ yet the 
difference would have been such as to produce no difference 
in the final result, as applicable to any physical investiga- 
tion in which the undetermined function is to be adapted 
to given physical circumstances. A simple instance will 
shew this. 

Let a (a?, y) = x. ^. Then we obtain the two following 
series for the integral: 

First, 
z ^ (f>{ax +y) +T2-2{6aV/ +4aVy +aV}. 

Second, 

The second value, it is easily seen, may be put under 
the form 

'f {ax +y) + —^ {GaVy* +4aVy +aV} 
-7^ y +4aa^ -V^a^a?]/^ +4aVy +aV} 
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= f {aoi +y) - j2^, {ax +y)* 



+ j^, {eo'ajy +4aVy +aV}. 



Now — r^-i (<KC +y)* is a iiinction of {ax +y) ; and, when 

we proceed to determine the form of the function '>fr{ax +y), 
which must be adopted, in order, with the other terms of 
the expression, to satisfy any condition prescribed hj ' 
physical considerations, we must at the same time take 

intoaccount— — — J {ox+yY as a function of {ax+y). We 

may therefore at once unite it with '^ [ax +y). Call their 
sum X (^^ +y)« Then the second solution is 

« = X («a; +3^) + j2^2 (6aV/ +4aVy +aV); 

precisely the same as the first, except that the symbol x 
stands in the place of <^. But in both cases the form of 
the function <f) or ^ will be determined by the consideration 
of satisfying the same physical condition. And therefore, 
whether we write ^ or ;^, we shall afterwards arrive at the 
same expressions; and therefore the two solutions, first 
and second, are in use identical. 

28. In many cases, when the form of solution has 
been obtained by the process of Article 25, the details will 

C2 
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be obtained most readily hj assuming the form of solution 
with indeterminate coelEeients. 



29. It is not intended here to go into any discussion 
of the more complicated forms of Partial Differential Equa- 
tions of the first order. The following however may be 
mentioned as flowing immediately jfrom the solutions which 
we have found. 

and the equation is immediately reduced to the form 
already treated. 

If the first side is e(aj)x^ ^^^^iyj^^-f » 

and the equation is reduced to the same form. 

In reference to physical investigations, the theory of 
Partial Differential Equations of the first order is princi- 
pally valuable as introductory to that of the second order. 
In this view, it is hoped that the instances here given are 
sufficient. 
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TREATMENT OP THE SIMPLEST PARTIAL DIFFERENTIAL 
EQUATION OF THE SECOND ORDER. 

30. The simplest Partial Diflferential Equation of the 
second order is 

— = 
dx . dy 

Integrating with respect to a?, and remarking that, by 
virtue of the reasoning in Article 19, there must be added 
to the integral an undetermined function of y, 

Integrating this with respect to y, and remarking that, 
for the same reasons, an undetermined function of x must 
be added, 



GEOMETRICAL INTERPRETATION OF THE SOLUTION. 



31. This equation and its solution admit of easy 
geometrical illustration. Keferring to Article 11, it will be 

seen that -7- is the trigonometrical-tangent of the angle 

between y and the curve-tangent, in the curye formed by 
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the intersection of a curved-surface with a plane parallel to 

zy. And the equation ^ — -r- = ^' ^' ^~ ( X ) ^ ^> denotes 

dz 
that, when x is changed, -^ undergoes no change; or that 

ay 

the inclination in the plane zy will be the same for a 
point more or less advanced in the direction of x than for 
the point under consideration. The curve-tangent may be 
absolutely higher or lower in the direction of «, but it will 
preserve the same inclination to y. As this applies to 
every point of the intersection-curve, it follows that the 
intersection-curve more or less advanced in the direction of 
X will be identically the same as that at the point under 
consideration; or that all sections by planes parallel to zy 
give the same curve, though perhaps at different elevations. 
And from this it follows (considering that, when two similar 
curves are separated in the direction of x, the slope from 
every point of one to the corresponding point of the other, 
in the direction of x, must be the same) that all sections by 
planes parallel to zx give the same curve, but not necessa- 
rily the same as those given by planes parallel to zy. Then 
the elevation of any point z may be considered as composed 
of these two parts; first, the elevation of the corresponding 
point of that zy curve whose plane passes through the 
origin, which elevation (since the form of the curve is the 
same for all values of x) may be called ^ (y) ; secondly, the 
elevation of the point z above the point last considered, 
which elevation (since the points are connected by a curve 
whose form is independent of y) may be called -^ (a?). 
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The whole elevation of the point z will therefore be 

^ (y) H-'^ (^)- 

32« These Amotions may be discontinnous, for the 
reasons stated in Article 22. 

33. In figure ^is a representation of one of the solids 
in which ^ {jj) and '^ (a?) are both continuous functions. 
It is scarcely necessary to remark that the forms of such 
solids may be infinitely varied ; thus (taking the vertex as 

origin of coordinates), if « = — V^ , or both sections be 
parabolas with different parameters, the curve-surface is a 
paraboloid with elliptic base; if « = — = - , the curve- 
surface is the paraboloid of revolution (firom which, figure 

2 was drawn) ; if « = '^ , the curve-surface is like a 

a 

limited portion of the interior of an annulus, or like a 

mountain-pass, &a &c. &c. 

In figure 3 is a representation of a curve-surface in 
which one function is continuous (the sections being similar 
parabolas), the other is discontinuous (each section being 
two sides of a triangle, the triangle being the same 
throughout). 

In figure 4 is a representation of a surface in which 
both functions are discontinuous (each section being two 
sides of a triangle, but the triangles in the plane xz being 
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different from those in the plane yz). It forms, in fact, a 
pyramid with trapezoidal base; two vertices of the trape- 
zium having the same value of cc, and two having the 
same value of y, and the vertex of the pyramid being above 
the intersection of the diameters of the trapezium. 

34. The solution of the equation -^ — -^ = a (a?, y) 
obviously is, 

z = <l>{y) +ylr (a?) +1 . . a {x,y). 

J y J X 



TREATMENT OF ANOTHER PARTIAL DIFFERENTIAL EQUA- 
TION OP THE SECOND ORDER. FIRST, BY CHANGE OF 
INDEPENDENT VARIABLES. 

35. To solve the equation j^^—^' j^ = ^ (^> y)* 

[This equation is the most important of all, especially in 
reference to those physical theories in which wave- 
transmission, of sound, of light, of water, &c., is explained 
by mechanical theory.] 

The best method of solving this equation is by the 
Change of Independent Variables. 

Let u = ax+y, v = ax—y 

which giTe a: = -^,y j. 

Consider 2? as a function of x and y because it is a 
function of u and v. 
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Then 

dz dz du dz dv dz dz 

T = :r^T" '^T^T = T^^ +:7-xa; 
dx du ax dv OM du dv 

d^z ^ d ldz\ __ d (dz\ du d ldz\ do 
da? ~" dx \dx) "" du \dx) dx dv KjoUc) dx 

d (dz dz \ . d fdz . dz \ 

= -v-(-y-xa -^--y-y^a xa +3- 13-xa + 3-xa xa 

du \au dv J dv \du dv / 

- «^ 4.2a^-^-?— +a'— • 
"" du^ du . fl?v dv^ ' 

dz ^ dz du dz dv _ dz 4. ^y_-|. 
dy "~ e?M dy dv dy du dv * 

d?z _ d (dz\ __ d (dz\ du d (dz\ dv 
dy^ " dy \dy) "" du \dy) dy dv \dyj dy 

d [dz ^ dz S\ . , ^ f^^ ^ dz ,\ , 
du\du dv J dv\du dv J 

d^z d^z d^z 



Hence 



du^ du . dv dv^ ' 



d?z 



da? dy^ " du . dv' 



And the original equation, divided by 4a', becomes 



_cPz 
du 



?z ^ 1 fu +v u — vN 

Tdv " 4^*^^ v"^^ ' "2~; • 
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Whence, by Article 34, 



or 



u = <l>{ax-y)+ir{ax+y) + ^,f^. j^. a (^, V)' 



where, after the integration, for u and v are tobeputaa5 +y 
and ax — y. 

36. When one of the independent variables is time, 
it will be seen by the same reasoning as in Article 24 
that the two undetermined functions represent two waves 
travelling in opposite directions. 



SECOND TREATMENT OF THE SAME EQUATION, BY SEPARA- 
TION OF THE SYMBOLS OF OPERATION FROM THOSE 
OF QUANTITY. 

37. A second method of solving the equation 

is founded upon the "separation of the symbols of opera- 
tion from those of quantity." This theory is, in fact, 
merely a convenient form for exhibiting the indubitable 
results of legitimate algebra; but it sometimes serves to 
suggest new methods of treating equations, which, when 
verified, are useful. 
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38. The equation 

may be written, almost without departure from ordinary 
notation, 



(s+«i)'-(^-^>' 



the connexion of the left-hand bracket with z which follows 
it being, however, not by multiplication but by differen- 
tiation. And, if 

we may, on the same principle, write it 
and therefore 

still on the understanding that the connexions on the left 
hand are not by multiplication but by differentiation. 

39. If we treated all the left-hand symbols as signs 
not of operation but of quantity, their product would be 
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or, pursuing the fanciful idea still further, 



d^ "^ d^*' 



and this suggests the idea that, if we follow up the true 
diflferential operations instead of the fanciful algebraical 
multiplications, we shall arrive at that result in its true 
diflferential meaning. And so, in fact, it proves. For, by 
actual differentiation, 

d fdz dz\ ^ d^z d*z 

dx \dx dy) ~" da? dx.dy^ 

d fdz dz\ __ d^z ^ %d^z ^ 

dy \dx dy) dx . dy dif ' 

and, adding the two lines, 

( d d\ /dz dz\ __ d'z j d^z 
\dx dy) \dx dy) da? djf^ ' 

Here the analogy with the algebraical operation has led to 
a real gain of convenience, by shewing that we can here 
break up the operation for a Partial Differential Equation 
of the second order into two operations for the first order. 

40. Thus, to solve the equation 

d^z ^d^z . . 
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we have first to solve the equation 

from which (X, Y) will be found; and then to solve the 
equation 

from which z will be found. 

Both operations are effected by the process of Article 25; 
the most convenient assumption however being 

u = ax +y, V = aa? — y. 
We shall return hereafter to the " separation of symbols." 



OTHER PARTIAL DIFFERENTIAL EQUATIONS OP THE 

SECOND ORDER. 

41. To solve the equation 

cPz , d*z 

or 

It does not appear possible to solve this equation gene- 
rally, except by the use of imaginary symbols. If we take 
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the second form of the equation (as written above), and 
apply the solution of Articles 35, &c. 

If, as is usual, we suppose <^ and '^ to represent real 
functions, it will be impossible to destroy the imaginary 
terms in this expression except by making -^^ the same as 
^, and thus giving up part of the generality of the solution. 
But the introduction of two undetermined functions may be 
thus obtained. Let 

^ = X+V— 1.0), -^rrrj^ —V —1 . fi), 

where x ^^^ ® ^® real functions. 
Then 

+n/ — l.fi)(y +a V — 1 .oj) —V — l.a)(y — aV -l.oj). 

Here the imaginary terms are destroyed, and there are 
two undetermined functions. 

It will be seen hereafter that, in some practical appli- 
cations of this equation under limiting conditions, the 
solution may be obtained more easily by specific process 
than by using the general solution. 

42. The following equation, which occurred to the 
author in the investigation of the movement of spherical 
waves of air (gravity being neglected), deserves attention 
for the nature of the substitutions made and the form of 
the result. 
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The equation given by the physical considerations is, 

d^z ^ t d /2z ,dz\ 
dt " ' dx\x dxj* 

The last term may be made simpler by assuming 

or y =s xj z. 

Differentiating this twice with respect to x^ 

d^y _ 9 , dz ^ 
da? "" dx^ 

A 1 rfV _ 20 dz ^ 

x' da? "~ X dx^ % 

whence the second side of the equation 



dx\x da?) 



For the first side, or -^a , we have 



d?z 
df 



_ d^ d^ (y\ _ _^ ^ (y\ 

■" drdx'Kx) " dx' de\x)' 



which, since ^ can be affected by -^ only because y is a 
function of ^, becomes 

d^n rfVV 

dxWdeJ' 
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and the original equation is changed to 

dx\x' df J ' dx\x' da?) 
Integrating with respect to x, 

1 ^ _a^ 4V _ -vrt 
x' de x'da? ~ ^ ^'' 

Solving this by the method of Article 35, 

y — a; . X («) +^ {at +ic) ^--^ [at —x) ; 
whence 



--^'(a«-aj). 



43. We believe that scarcely any other equations than 
these, occurring in physical investigations, have been solved 
in a finite form. Several equations have been solved in 
the unsatisfactory form of infinite series, of which the con- 
vergence is not always assured ; but it is not the object of 
the present Treatise to enter on them. 

44 The following equations (probably among many 
others occurring in physical inquiries), which have pre- 
sented themselves in the author's investigations, have not 
been solved. 
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In investigating the vertical motion of a horizontal 
wave of air, on the supposition that the elasticity is propor- 
tional to the density, and that gravity is constant, we 
arrive at an equation of this form : 

d^X dX 7 2 ^ _ ^ 

which we have not succeeded in integrating generally in a 
finite form. 

In investigating the radial motion of a spherical wave 
of air whose center is the center of the earth, supposing the 
elasticity to be as the %*** power of the density (where n 
must be greater than a certain quantity which is > 1, in 
order to make the whole mass of atmosphere finite), the 
resulting equation takes this form : 

|^_(n-l)aV /d5 ^^ _n^fdE R. _^^^^ 
[ r )dr\dr rj r \dr rj dtr 

This equation may be much simplified, but in any 
form we have not rendered it integrable. 



FURTHER CONSIDERATION OF THE SEPARATION OF THE 
SYMBOLS OF OPERATION FROM THOSE OF QUANTITY. 

45. In Article 39 attention was called to the method 
of solving the equation then under consideration by the 
"Separation of the symbols of operation from those of 
quantity." This principle may be applied to an extensive 

D. £. D 
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system of equations, which are however, at present, matters 
rather of curiosity than of physical value. It may be 
sufficient here to indicate one class. 

If the equation 

can be expressed as 
[L -"^)^(^ -5^)x&c.ton terms x^ = a (x, y), 

then the process of Article 40 can be applied with so little 
difficulty that it appears unnecessary for us to delay 
further on it. In the successive changes of Independent 
Variable, it will be found convenient to take the factors in 
successive pairs; thus, for the effect of the first pair of 
factors, where 

(s-«|)(s-»|)(^•'^-«(^.^)• 

assume ax ■\-y = w, hx +y = v, 

and proceed as in Article 25. 

46. There is one exceptional case, however, well 
brought to notice by this method of treatment, which 
merits further attention ; namely, that in which two factors 
a, b, are equal. To take the simplest case, suppose 

^ _2 ^^ «^^ - 
dx^ dxdy dy^ ' 
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or 



\dx dy)\dx dy ) 



Here we may assume ax +y = w, ea; -^-fy = v, e andy* 
being any constants whatever whose proportion is not a: 1; 
the value being admissible for either. 



Then 



d(X,Y) d(X,Y) d{X,Y) 
dx du dv 



dy du dv "^ ' 

and the equation 

becomes 

whence (X, F) = ^ (w). 

And, putting for (X, Y) its value, 



(|-"|)" = '^(")5 



or, by the same substitution which was made in regard to 
the differential coeflBicients of (X, F), 

D2 
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integrating this, and remarking that ^(u) possesses the 
properties of a constant as regards integration with respect 
to r, 

or 



z 



= [r=^fy +7^/^^ ) • * (y +'^) +^ (y +«*)• 



Observing that a constant factor may be supposed to be 
included in the form of an undetermined function, we may 
express z in either of the forms, 

{gy +hx)x<l> (y +ax) ^--^ (y +ax) ; 

^x0« (y +^) +'^« (y +^) ; 

all which, as applicable to the satisfaction of any assigned 
conditions, will be found to be identical, on employing the 
reasoning of Article 27. 

47. In pursuing this subject, the signs / and d may be 
considered as reciprocal, or J^cfK In some cases, equa- 
tions may be valued by carrying symbols, such as 

dx dy^ 

^hich hold the position of inultiplier on one side of the 
equation, to the position of divisor on the other side, or 
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aa multiplier in the form 

d 



i^^a -1 
\dx dyl 



And in this shape they may be substituted in abstruse 
and general theorems. For instance, to take a case which 
(in this theory) is almost elementary, we know that 

If n = — 1, this becomes, reversing the sides of the equation, 



e 



or 



(i+«) '"=*'"//*'"• 



Bat the equation 



dz dz _ 

dx dy ' 

d d 



(d , d\ 

may, in accordance with these principles, be put in the 
form 



\dx dy) ' 
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and the second side agrees with the first side above, 
provided that e -7- be put for a and treated as a constant. 

Then we shall have 

<2 



J X 



a solution to which, in some cases, a real meaning maj be 
given. 

48. This principle, as a purely algebraical and 
symbolical process, possesses very great power, and leads to 
very remarkable results. But the reader cannot fail to 
observe that it carries with it no evidence whatever for 
the validity of results, (such as is conveyed by the opera- 
tions of quantitative algebra, or by the steps [properly 
pursued] of the difierential calculus), for which it must 
rely on subsequent verifications. As aiding the appli- 
cation of Partial Differential Equations to physical 
investigations it possesses little value. The further 
examination of it would therefore be out of place in 
this Treatise. The student who desires to follow it up 
will find much information in Boole's Treatise on Drffk" 
rential Equations, Gregory's Examples, and similar works. 
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TREATMENT OP PARTIAL DIFFERENTIAL EQUATIONS WHEN 
THE SOLUTION IS LIMITED BY PREARRANGED CON- 
DITIONS. 

49. We have seen in Article 44 that there are some 
Partial Differential Equations of which the solutions 
cannot be obtained by known methods, so long as the 
utmost generality is required in the solutions. But in 
some cases, by attaching limiting conditions to the solu- 
tions, the same equations may be integrated with ease ; as 
will be seen in the following instances. 

60. In considering the tidal disturbance of water in an 
equatoreal canal round the earth, retarded by friction pro- 
portional to the velocity, we arrive at the following equation 
{x measured along the canal, X the horizontal disturbance 
of a particle) : 

= Hsin {it —mx) —f-jr +a 



df ^ ' -^ dt dix?' 

where JT, f^ a*, are constants depending on the moon's 
attraction, the coefficient of friction, and the depth of the 
sea. This equation cannot be solved generally. But, 
remarking that the only part of the solution which has any 
interest for us is that which follows the same law of pe- 
riodicity as the lunar motions, we may assume 

X = u4 sin {j,t —mx) +£cos {it — wa?), 
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where A and B are constants to be determined. On 
substituting this, the values of A and B and the expression 
for X are found without difficulty. 

51. Before adverting to a specific solution of the 

equation ^ +«* T-a = 0, we will make the following 

remark on the general solution. When an exponential 
form is adopted for the functions in the general solution, 
Article 41, the general form may be retained conveniently: 
but in other cases it is usually necessary to expand tha 
functions. Performing this process according to the ordi- 
nary algebraical rules, the general solution becomes 

^ («) -^" {x) ^' +f' " {X) Y^^ -&c 

+ t'(a')f-t"'(»)i^ + &c 

which, as will be found on trial, satisfies the equation. 
A solution by series, however, can rarely be considered as 
quite sufficient. 

52. In considering the motion of ordinary small waves 
in a sea of uniform depth {x horizontal in the direction in 
which the wave is going, y vertical measured upwards from 
the bottom), we find 

^ ^= 
dy* do? "" 

The cumbrous form of the general solution of this 
equation wpuld make it almost useless to attempt to apply 
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it -to the special case. But the nature of the case permits 
us to assume 

X = Iicos{nt '~mx)+S sin (nt —mx), 

m 

where JS and 8 are functions of y. Substituting, we have 

whence 

B = Ce'^+De^, 8 ^ C e"^" +1)' €^ ; 

and, after some reductions peculiar to the problem, 

X ^A (e*"^ +€"^«') cos {nt -mx -B). 



FINAL DETERMINATION OF THE FUNCTIONS WHICH ARE 
UNDETERMINED IN THE GENERAL SOLUTION. 

53. In nearly the whole of these solutions, our r©^ 
suiting expression for the quantity which it is our object 
to find is accompanied with undetermined functions. It is 
now an object of great importance to determine those func- 
tions. Although in reality the process for doing this may 
be stated in brief rules, yet it will be better understood 
from an exhibition of its application in one or two actual 
instances than in any other way. 

54. Example 1. In the consideration of the tidal 
motions of water without friction in an equatoreal channel 
round the earth, X being the horizontal displacement at 
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any time ^ of a particle whose ardinate measured round the 
channel &om. a fixed origin was x, the equation is 

where it is double the moon's hour-angle from the origin, 
and mx is double the longitude-angle of the point x from 
the same origin ; and where c* is a constant depending 
on the depth of the water. This equation can be com- 
pletely solved : its result is 

^ — 2 «_ ^ -sin {it -^mx) +(f> {ct -\-x) 4--^ {ct —x) ; 

which we shall write 

C sin {it —mx) +(f> {ct +x) ^-^ {ct —x). 

Each of the three terms represents a wave, but they 
are waves of different characters. The first wave corre- 
sponds in its velocity with the lunar force Hain {it — mx) 
which produces it: for this the author has introduced the 
term forced wave; it depends entirely on the lunar force, 
and would not exist if there were no lunar force. The 
second wave is travelling in the direction opposite to 
the moon's apparent diurnal motion, and the third wave in 
the same direction as the moon's apparent diurnal motion j 
for these the author has introduced the term free waves; 
their velocity is independent of the moon's motion; the 
two waves, or either of them, may or may not exist; they 
are entirely independent of the moon's action, and the 
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q&estion of their existence is generally independent of 
the moon*s action and of the forced wave. This is the state 
of things when I^HER is no material limit to the length of 
the canal. 

55. But suppose that the canal is limited hj the ob- 
stacle of land at its two ends (as in the Mediterranean Sea, 
all minor gulfs, &c., being put out of consideration). Then 
the forced wave alone is not sufficient for a solution. For, 
the existence of the terminal obstacles (whose ordinates we 
will call a and h) requires that, at those obstacles, the hori- 
zontal displacement of the water be nothing at all times; 
and therefore when a? = a or a? =s J, X must = whatever 
be the value of t 

Now this cannot hold with the forced wave alone; for, 
the values which the first term assumes for a? = a and a? = J^ 
namely, Cain {it— ma) and C ain {it —mb), will not vanish 
for all values of t. We must therefore have recourse to the 
two undetermined functions; or, in other words, we must 
now necessarily introduce the two free waves as an indis- 
pensable part of the solution : and we must determine their 
elements so that the two conditions, of X always = when 
x=sa and when x ^b, shall be satisfied. 

56. The term expressing the forced wave is a simple 
periodical term having it in its argument. Therefore, in 
order to destroy this at all times for definite values of aj, 
the free waves must be expressed by simple periodical 
terms having: it in their arguments. Therefore the two 
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undetermined terms ^ {ct +ixt) and '^ {ct —a?) must have 
the forms 



7 



and 



A^ sin -j- (c< + a;) f +B^ cos -^- (c< +a;) \ , 



-4, sin j- {ct —x) y +5, cos ] - {ct —x) \ ; 



and the entire expression for X will have the form 

CBm{{t'-mx) 

• • 
4-ili sin (tif +-0?) + B^ cos (tV +- a?) 

• -• 
+-4, sin (i^ — x) +5j cos («^ — x)-, 

or, expanding the trigonometrical terms, and putting e for 
A^ +^„/ for ^, - J„ g for 5^ +5,, A for 5, -J5,, 



sinti^xj 



to; 



C cos 9wa? +6 cos hA sin 

c 



. %x\ 



t<xl —(7 si] 



«o; tXl 

+ COS t<x ^ — (y sin ?waj +/ sin Vg cos — h . 



Making this = when x=^ a and when x = J, and 
writing a for — and /8 for — , we have the four equations, 

= (7cos?wa +ecosa +Asina, 
=5 — C7sin wa +/sin a +5^ cos 0, 
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=» (7cos wJ +ecos/8 +Asin/8, ^ 

= —C sin mb +/siii fi +g cos 13, 

Solving these equations, 

e sin (/8 —a) = (7(— cos ma sin ^8 +cos mb sin a), 

/sin (^ — a) = (7 (— sin ?na cos /8 +sin wJ cos a), 

^ sin (^ — a) = C{ sin wa sin ^ —sin wJ sin a), 

A sin Q8 — a) = C7 (cos wa cos ^ —cos wJ cos a). 

Substituting, and restoring the original notation,. 

'''^ * (c^m»-i«)8inr-(J-a)l 

sin (ii —ma) sin |- (J —a;) > +8in(t^ — wJ)sin \-(ps — a)[ . 

Thus the terms, undetermined in the general solution, 
are now completely determined so as to satisfy the special 
conditions of the problem. The entire expression for X, it 
will be seen, satisfies the original partial difierential equa- 
tion, and also makes X = at all times when X ~ a and 
a; = J. 

57. Example 2. In the orditary problem of vibrating 
musical strings, where no force is supposed to act after the 
string has been put in motion, x being measured longi- 
tudiu,ally along the string and z being the smaU transversal 
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dUplacement of anj point, the equation for z is 

d^z , d^z 

where a =gL, L being that length of the same string whose 
weight will represent the tension of the string: and its solu- 
tion is 

z = ^ {at +x) +'^ {at — a?). 

Here the solution is expressed entirely hj the two 
undetermined functions. If the initial circumstances of 
disturbance (namely, the displacement and the velocity of 
every point) were given, those functions could be deter- 
mined, as will be shewn below. But even without an 
absolute knowledge of the forms of the functions, some 
very important properties of the solution may be ascertained 
thus. 

58. Let the ordinate of the near end of the string be 
0, and that of the further end L Then, whatever be the 
yalue of <, « is when a? = or when x = I That is, 

{at) ^^Ir {at) = 0, 
<f>{at+T) +'^ {at -I) = 0. 

The first equation shews that yfr must always =s — ^. 
Thus one of the undetermined functions is already elimi- 
nated (the other will be eliminated from consideration of 
the initial circumstances) r 

Substituting in the second equation, 

<l>(at+l) -'<f>{at-l) = 0. 

Now at may, with changes of time, receive any value 
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whatever. And if a^ — Z = j (which makes at +1 = q +2l), 
q may have any value whatever. And thus, q having any 
value whatever, we always have the equation 

<i>{q+2l) =^ (f>{q); 

that is, the form of the function ^ must be such that, on 
increasing or diminishing the quantity under the bracket 
by 2Z, the same value of the function will be reproduced. 
Now in either term of the expression for the displacement 
of the particle at x, namely 

ff} [at+x) — ^ {at —x), 

the increase by 2l of the quantity under the bracket may 

2l 
be effected by increasing at by 2Z, or by increasing ^ by — . 

a 

2l 
Therefore, every time that we increase ^ by — , we find 

the particle x with the same displacement as before. And 
this applies to every particle of the string. Consequently 
the string performs a complete vibration, returning again 

... .21 

to the same position, in the time — . 

It is true that, as will shortly be seen, under certain 
circumstances the string may perform two or more com- 
plete vibrations in the same time, but it always returns to 

2l 
the same state after the interval — : and this is thQ only 

a "* 

interval of which it can be asserted that a return to the 
same state necessarily occurs in that time. 
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The reader may exercise his ingenuity in proving that 
the motion of the string consists in the movement of two 
waves, each of which travels from end to end with velocity 
a, and is reflected at each end. 

59. Putting, for the moment, for the quantity 
affected by the functional symbol 0, it appears that ^ {&) 
is a finite periodical function, going through its changes 
and returning to the same value when is changed by 2?. 
A trigonometric function of sines and cosines of an angle 
goes through all its changes and returns to the same value 
when the angle is increased by 27r. Hence it appears that 
^ (G) is similar in its general character to a trigonometric 

function of sines and cosines of — —j— , ov j 0. Any 

finite periodical magnitude may be represented with any 
assigned degree of approximation by a series of integral 

powers and combinations of the sine and cosine of j 5 ; and 

these powers may be converted into sines, and cosines of 
multiple arcs. Thus the function (f> {6) may be represented 

by ^^ 

-4^ sin J ^ +A^ ain-j-O + &c. 4--4» sin -y ^ +&c, 

+B^ cos J 6 +5, cos -J- 6 +&C +5, cos -y +&c.; 
which may be conveniently expressed by 

S(A8in^e) +2(Acos^5); 
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and thus the value of z will be 

2 \A^ sm — ^ *-> +S JAi cos — ^—. ^ • 

-S lA^ Sin — ^—^ '-)' -2 j-Bn COS — ^ ^-V , 



or 



^ f ^ A mrat . mr \ 

2 ( 2-4^ cos — y- Sin -J a; 1 



—2 ( 2-Dh sm —J- sm -y- a? I ; 

dz 
and the value of -7- will be 

^ f2mra . . W7ra< . nir \ 
—2 ( — ^ — -4» sin — J— sin -y- x\ 

^ /2w7ra T^ WTraf . nir \ 
—2 ( — J- B^coa —J sm -=- wj . 

60. Suppose now that the initial circumstances of 
displacement (Z) and velocity {Z') are given for every 
point of the string, and that from these we desire to deter- 
mine all the subsequent motions* We must now make 
^ = 0, and we have 

Z = X (^A^ sin-y- a; j , Z' = —2 (-y- n Bn sin -t-x\ . 

D. £. £ 
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And our immediate object is, knowing the values of 
Z and Z* for every value of a?, to determine the values of 
u4j, ^j, &c., Byy J5j, &c. This will be done by means of 

the following general formula. If we multiply sin -^x 

by sin -y- a?, where m is any integer different from w, and 
integrate from a: = to a? = Z, we have for product 



-cos ^ X --cob'—^-I-x, 



for general integral 



I . (m —n)ir I . (m +n)7r 

sin -i =—L- X ——7 ; — r— sin-^ — -. a?, 



2(m — 7i)7r Z 2 (m +w)7r Z 

and for definite integral 0. 

But if w = w, we have for product 

1 1 2mr 
- — cos — 5— a?, 

2 2 1 ' 



for general integral 



X I . 2/i7r 



and for definite integral -. 
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Hence, if we put 8 for the integral with respect to x 
taken from a? = to a; = Z, 



/Sf^ZsinyO?) = A, I; 
8(Zain-j-xj — A^l, &c. 
S (Z' sin -j-ic] = — 7rawJ5j; 
8 (Z* sin — T— xj = --wan J5„ &c. 



61. In a given instance, in which the initial displace- 
ments and velocities are given not by formula but by 
numerical values, the following process will apply. Sup- 
pose we limit the multiples of -j 0j where sines and cosines 
are to be used, to the first six, namely, 



TT^ 27r>, Sir ^ ^ir^ hir ^ 6w ^ 
T^^ T^' T^' T^' T^' T^' 



which will suffice in almost any conceivable case. Divide 
the length I into 120 equal parts, and suppose that Z and 
Z are known numerically for the middle of each of these 
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nir 



parts. For the first of these middles, -^ a? = 45'; for the 

second it is 3x45'; for the third it is 5x45'; and so on. 
And for each of these parts, the integral is, in fact, taken 

through the extent -— . Hence, the formula above gives 

the following numerical rule : — 

Multiply the successive values of Z by sin 45', sin 3x45', 
sin 5x45', &c.: the sum of products = 120 A^. 

Multiply the successive values of Z by sin 2x45', 
sin 6x45', sin 10x45', &c.: the sum of products = 120-4^. 

Multiply the successive values of Z by sin 3x45', 
sin 9x45', sin 15x45', &c. : the sum of products = 120-^3. 

And so to Aft. 

And a similar process (with proper changes) determines 
J5j, J5j, &c., from the given values of Z*. 

Then substituting these in the formula at the end of 
Article 59, the circumstances of all subsequent motion are 
completely obtained. 



62* If A^y or JSj, or both, are found to have real 
values, and -4^, A^y &c., J5,, B^y &c., have no real values, 
tlien the vibration is that of a simple line of sines, de- 
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pending, for time, on the angle — p . If only A^ and jB^, 
one or both, have real values, then the vibration is that of 
two lines of sines ; separated at the point where -y- a? = tt, 

or where a? = -5 , thus forming a permanent node in the 

middle of the length ; and depending for time on the angle 

— J— , and thus vibrating in half the time of the former 

vibration. Similarly if A^ and B^ only have real values, the 
string is divided into three equal parts by two nodes, and 
the vibrations are made in one-third of the time of the 
fundamental vibrations. 

The connexion, of these results, with the results of the 
experience of the senses as to the musical tones produced, is 
one of the most important points in the Acoustical Theory 
of Music. 



63. The two instances which we have given will 
probably suffice, better than any rules, to shew the use 
that is to be made of the undetermined functions in the 
solution of a Partial Differential Equation. 
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CONSIDERATIONS ON THE NECESSITY OF CONNEXION 
BETWEEN THE NUMBER OP UNDETERMINED CONSTANTS 
AND THE ORDER OP THE EQUATION (iN SIMPLE DIF- 
FERENTIAL equations), and between the NUMBER 
OP UNDETERMINED FUNCTIONS AND THE ORDER OP 
THE EQUATION (iN PARTIAL DIFFERENTIAL EQUA- 
TIONS). 

64. In solving a simple differential equation of the 
first order, we usually arrive at a solution containing one 
undetermined constant. But it is not always so. There is 
a remarkable class of solutions called "particular solutions," 
which not only have not, but which cannot have, an un- 
determined constant. Thus, if we express in algebraical 
language the problem "To find the equation to the curve 
whose tangents possess this property, that a perpendicular 
from a given point upon the tangent has a given value ;" 
we have, for general solution, a straight line with one 
undetermined constant; and for particular solution, a circle, 
which from the nature of the case does not admit an un- 
determined constant The connexion, therefore, between 
the order of the equation and the number of undetermined 
constants in its solution is not invariable. 

65. But if we assume a formula as solution to a dif- 
ferential equation which we propose to form, the said 
formula containing symbolical constants, we can pre- 
determine that any one of those constants (if the equation 
be of the first order), or any two of those constants (if the 
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equation be of the second order), shall be eliminated from 
the differential equation to be formed, and shall therefore 
be undetermined when that equation is solved. Thus, to 
take a very simple instance, let y = ax +i be assumed to 
be the solution of a differcDtial equation, with the constant 

a undetermined. Differentiating it, -^ = a ; substituting 

this in the given equation, y— a?-^— J = 0, the dif- 
ferential equation required. This process evidently can be 
extended to equations of any order. 

66. It appears thus that the existence of a soluble 
differential equation does not necessarily imply an unde- 
termined constant in the solution; but the existence or 
assumption of an undetermined constant in a solution can 
always be represented by proper form of a differential 
equation : and that this theorem applies to equations of 
jevery order. 

67. If now we apply similar considerations to Partial 
Differential Equations, we are led to the following con- 
clusions. 

68. First, on examining the solution in Article 18, 
&c., it will be seen that, in the cases in which the solution 
of a Partial Differential Equation is really, effected, it has 
been done by reducing it to an ordinary integration, just 
as would be done for a Simple Differential Equation. This 
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seems to indicate a probability that there may be cases in 
which the solution of a Partial Differential Equation has 
limitations analogous to those of a Simple Differential 
Equation, and therefore has not, and cannot have, an un- 
determined ftinction in its solution. Upon this, however, 
there is not at present any distinct evidence, and the mat- 
ter is indicated as one which appears to merit exami- 
nation. 

69. Second, if we assume a formula containing one 
function ^(w?), where t^ is a definite fimction of x and y, 
and where <^ is intended to be undetermined, we can always 
produce a Partial Differential Equation of the first order, 
of which the assumed formula will be the solution. 

For we have 

Zf which contains (f>{w) ; 

and we can form 

dz 

-,- , which will contain (}>{w) and (}>'{w) ; 



•5- , which will contain ^{lo) and ^'{lo) ; 

and from these three equations we can eliminate fl>(w) and 
<l> {w), leaving one equation between z, x, y, -7- , and 



dx 



dz 

-^ , of the form required. 
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70. But, third, if we assume a formula containing two 
functions (f>{w) and '^(5), where w and s are definite func- 
tions of X and y, and where ^ and yfr are intended to be 
undetermined, we cannot always produce a Partial Dif- 
ferential Equation of the second order, of which the as- 
sumed formula will be the solution. For, we have 

z, which contains <l){w) and ylr(s) ; 
and we can form 

dz 

-r- , which will contain <^(tt?), ^'{w)y -^(a), •^'(*) ; 

dz . . 

-T- , which will contain ^(t^), ^'(^)j "^Wj "^'W > 

d'z 

-^ , which will contain <^(w>), ^'(*^)j ^"(^)> "^W* 



docdy 
d'z 



, which will contain the same six functions ; 



dy' 



, which will contain the same six functions. 



Here we have six equations, and no more, with six 
functions on the other side of the equations. We cannot 
therefore eliminate these six functions, so as to leave an 
equation between 

dz dz d^z d^z d^z 

D. E. F 
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And thtw we cannot always produce the Partial Dif- 
ferential Equation required. The same remark, it will be 
found, applie. more Btrongly to eqtiations of a higher 
order. 

71. It seems not impossible that the failure of at- 
tempts to solve the equations in Article 44 may have some 
connexion with this inability to establish a perfect con- 
nexion between the order of the equation and the number 
of undetermined functions. The subject appears to be 
worthy of greater attention than it has received. 
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